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I. INTRODUCTION 

Can thermal fluctuations be rectified? Ever since Maxwell Q raised this question with his famous thought experi- 
ment involving a so-called Maxwell demon, it has been the object of debate in both the thermodynamics and statistical 
physics community. The mainstream opinion is that rectification is impossible in a system at equilibrium. Indeed the 
property of detailed balance, which was discovered by Onsager 0, and which turns out to be a basic characteristic of 
the steady state distribution in area preserving time-reversible dynamical systems (in particular Hamiltonian systems) 
, states that any transition between two states (defined as regions in phase space of nonzero measure and even in 
the speed) occurs as frequently as the time-reversed transition. The separate issue, first introduced by Szilard jJl , of 
involving an "intelligent observer" that tracks the direction of these transitions, making possible the rectification by 
interventions at the right moment, has a contorted history of its own. It turns out that the engendered rectification 
is offset by the entropic cost of processing (and more precisely of erasing) the involved information Another more 
recent debate involving entangled quantum systems jQ] is still ongoing. 

Apart from the fundamental interest in the subject, a number of recent developments have put the issue of recti- 
fying thermal fluctuations back on the agenda. First, we mention the observation that thermal fluctuations can in 
principle be rectified if the system under consideration operates under nonequilibrium conditions. The past decade 
has witnessed a surge in the literature on the subject of the so-called Brownian motors 0- Such motors possibly 
explain, amongst other, phenomena such as transport and force generation in biological systems. Second, our ability 
to observe, manipulate or even fabricate objects on the nanoscale prompts us to look into new procedures to regulate 
such small systems, possibly by exploiting the effects of thermal fluctuations in a constructive way. 
Even though several constructions have been envisaged to discuss the issue of rectification in more detail, including for 
example the Smoluchowski-Feynman ratchet , the issue of thermal fiuctuations in a system with nonlinear friction 
Pi and the thermal diode no exactly solvable model has been put forward. In this paper, we will present two 

fully microscopic Hamiltonian models, in which the rectification of thermal fluctuations can be studied in analytic 
detail. Versions of the first model have appeared in the literature for some time under the name of Rayleigh piston 
or adiabatic piston [T^ ITsL IT^ . see also The second model, to which we will refer as thermal Brownian motor, 
was introduced in a recent paper |l6j |. Both models involve a small object simultaneously in contact with two infinite 
reservoirs of ideal gases, each separately at equilibrium but possibly at a different temperature and/or density. A 
Boltzmann-Master equation provides a microscopically exact starting point to study the motion of the object. As a 
result of the rectification of the nonequilibrium fluctuations, the object acquires a systematic average speed, which 
will be calculated exactly via a perturbative solution of the Boltzmann equation, with the ratio of the mass of the gas 
particles over that of the object as the small parameter. 

The organization of this paper is as follows. We start in section^]by reviewing the general framework and the type of 
construction for the Brownian motor that we have in mind. The main technical ingredients are closely related to the 
so-called l/f2-expansion of van Kampen In section UTTl we turn to a detailed presentation and discussion of the 
adiabatic piston. The rectification has in this case been investigated to lowest order by Gruber and Piasecki We 
present a streamlined derivation allowing to go two orders further in the expansion. Next, in section Hvl we discuss 
the more surprising thermal Brownian motor in which the motion derives from the spatial asymmetry of the object 
itself 16]. We again calculate the three first relevant terms in the expansion of the average speed. Finally, in section 
Ivl the obtained analytic results are compared with a direct numerical solution of the Boltzmann-Master equation and 
with previous molecular dynamics simulations |16|. 
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II. EXPANSION OF THE BOLTZMANN-MASTER EQUATION 



Consider a closed, convex and rigid object with a single degree of freedom, moving in a gas. To obtain a micro- 
scopically exact equation for the speed V of this object, we will consider the ideal gas limit in which: 

(1) the gas particles undergo instantaneous and perfectly elastic collisions with the object, 

(2) the mean free path of the particles is much larger than the linear dimensions of the object (regime of a large 
Knudsen number), 

(3) the (ideal) gas is initially at equilibrium, and hence at all times: the perturbation due to the collisions with the 
object are negligible in an infinitely large reservoir. 

With these assumptions, there are no precoUisional correlations between the speed of the object and those of the 
impinging gas particles, hence the Boltzmann ansatz of molecular chaos is exact .IZj. In fact, since the collisions with 
the gas particles occur at random and uncorrelated in time, the speed V of the object is a Markov process and its 
probability density obeys a Boltzmann-Master equation of the following form: 



H^(F|F') represents the transition probability per unit time to change the speed of the object from V' to V. Its 
detailed form can be easily obtained following arguments familiar from the kinetic theory of gases. 
To construct a model for a Brownian motor, two additional ingredients need to be introduced. First we have to 
operate under nonequilibrium conditions. This can most easily be achieved by considering that the object interacts 
not with a single but with two ideal gases, both at equilibrium in a separate reservoir, each at its own temperature 
and density. The physical separation (no particle exchange) between the gases can be achieved by using the object 
itself as a barrier (adiabatic piston) or by assuming that the object consists out of two rigidly linked (closed and 
convex) units, each moving in one of the separate reservoirs containing the gases. Second, we need to break the 
spatial symmetry. In the adiabatic piston this is achieved by the asymmetric distribution of the gases with respect to 
the piston. In the thermal Brownian motor, at least one of the constitutive units needs to be spatially asymmetric. 
With these modifications in mind, we can still conclude that Eq. ^ remains valid, but the transition probability is 
now a sum of the contributions representing the collisions with the particles of each gas. 

With the ingredients for a Brownian motor thus available, we expect that the object can rectify the fluctuating force 
resulting from the collisions with the gas particles. Hence it will develop a steady state average nonzero systematic 
speed, which we set out to calculate analytically. Unfortunately, an explicit exact solution of Eq. cannot be 
obtained even at the steady state, and a perturbative solution is required. Since we expect that the rectification 
disappears in the limit of a macroscopic object, a natural expansion parameter is the ratio of the mass m of the 
gas particle over the mass M of the object. More precisely, we will use e = y^m/M as the expansion parameter. 
In fact this type of expansion is very familiar for the equilibrium version of the adiabatic piston, namely the so- 
called Rayleigh particle. It has been developed with the primary aim of deriving exact Langevin equations from 
microscopic theory and culminated in the more general well-known l/f2 expansion of van Kampen [10|. With the 
aim of streamlining this procedure for the direct calculation of the average drift velocity, with special attention to 
higher order corrections, we briefly review the technical details. First it is advantageous to introduce the transition 
probability W{V';r) — W{V\V'), defined in terms of the jump amplitude r = V — V, since the latter jumps are 
anticipated to become small in the limit e — > 0. One can then rewrite the Master equation as follows: 



A Taylor expansion of the transition probability in the first integral of Eq. ((SJ with respect to the jump amplitude 
leads to an equivalent expression under the form of the Kramers-Moyal expansion: 




(1) 



dP{V, t) 

di 



W{V - r; r)P{V - r, t)dr - P{V, t) / W{V] -r)dr 



(2) 




(3) 



with the so-called "jump moments" given by 




(4) 



Since the change in the speed of our object of mass M, i.e., the jump amplitude r, will, upon colliding with a 
particle of mass m, be of order = m/M, the Kramers-Moyal expansion appears to provide the requested expansion 
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in our small parameter. However, the parameter M will also appear implicitly in the speed V . Indeed, we expect 
that the object will, in the stationary regime, exhibit thermal fluctuations at an effective temperature T^ff, i.e., 
^Miy"^) = \kBTcB_- To take this into account, we switch to a dimensionless variable x of order 1: 



The explicit value of T^s will be determined below by self-consistency, more precisely from the condition {x^) 1 to 
first order in e. The probability density P{x, t) for the new variable x thus obeys the following equation: 

71—1 ^ ^ 

with rescaled jump moments, An{x)^ defined as 



A^{x) = [j^Jl a^x). (7) 

Equivalently, and of more interest to us, the following set of coupled equations determine the moments {x") = 
/ x''P{x,t)dx: 

dt{x) = {A,{x)) 

dt{x^) = 2{xAi{x)) + {A2{x)) 

dtix"") = 3{x^Ai{x)) +3{xA2{x)) + {A3{x)) 

dt{x^) = 4(0:3^1 (x)) + 6{x^A2{x)) + 4{xA3{x)) + {A^ix)) 

dt{x^) = 5{x'^Ai{x)) + lQ{x^A2{x)) + lQ{x'^A3{x)) + 5{xAi{x)) + {A^ix)) 

dt{x^) = eix'^Aiix)) + 15{x^A2{x)) + 20{x^A3{x)) + 15{x^A4{x)) + 6(0:^5(0;)) + (^6(0;)) 

(8) 

The exact solution of this coupled set of equations is as hopeless and equally difficult as the full Boltzmann-Master 
equation. However, a Taylor expansion in e shows that the equations are no longer fully coupled and the calculation 
of a moment up to a finite order reduces to an in principle simple (but in practice tedious) algebraic problem. 



III. THE ADIABATIC PISTON 



A. Motivation 



In Fig. n we have represented in a schematic way the construction of the Rayleigh piston and of its nonequilibrium 
version known as the adiabatic piston. We concentrate here on an (infinite) two-dimensional system, for reasons of 
simplicity. The piston is considered to be a single "flat" particle of length L and mass M with a unique degree of 
freedom, namely its position x along the horizontal axis. Since the piston has no internal degrees of freedom, it can not 
transfer energy by "hidden" microscopic degrees of freedom. The absence of a corresponding heat exchange prompted 
the use of the name "adiabatic piston". The piston is moving inside an infinite rectangle separating the gases to 
its right and left from each other. These gases are initially taken separately in equilibrium, but not necessarily at 
equilibrium with each other. In the thermodynamic context of a macroscopic piston, this construction is an example 
of an indeterminate problem, i.e., the final position of the piston can not be predicted by the criterion of maximizing 
the total entropy, since it depends on the initial preparation of the gases ^| , see also . The case of interest to us 
is when the mass of the piston is not macroscopically large, i.e., finite e = y^m/M. When operating furthermore at 
equilibrium, this Rayleigh piston provides an exactly solvable model, allowing for example, the rigorous derivation of 
a linear Langevin equation appearing as the first nontrivial limit of the Boltzmann-Master equation in the limit e — > 0. 
When the left and right gases are not at equilibrium, but exert equal pressure on the piston, the model becomes an 
example of a Brownian motor, which is able to perform work by rectifying pressure fluctuations p^ . In doing so, the 
single degree of freedom x also plays the role of a microscopic thermal conductivity, an issue that is quite relevant 
to other models of Brownian motors 18]. Since this model is essentially one-dimensional and the related calculations 
are relatively simple, we include it in this paper to illustrate the calculation procedure and at the same time to derive 
novel results for the average drift speed up to order e^. 
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Figure 1: The adiabatic piston. 



B. Presentation of the model 



The ideal gases in the right and left compartments, separated from each other by the piston, are each at equilibrium 
with Maxwellian velocity distributions at temperatures T\ and T2, and with uniform particle densities p\ and p2, 
respectively. Since we are mainly interested in the rectification of fluctuations, we will focus on the case of mechanical 
equilibrium with equal pressure on both sides of the piston, i.e., piTi — p2T2- 

The motion of the piston is determined by the laws of Newton. Hence its velocity only changes, say from V to 
when it undergoes a collision with a gas particle, its (x-component of the) velocity going from to v^- Conservation 
of energy and momentum determines the post-coUisional speeds in terms of the pre-collisional ones: 



1 

—mv. 
2 



2 2 

mvrc + MV = mvL + MV' 



1 

-MV' 
2 



implying: 



V 



2m 



m + M 



(9) 



The transition probability then follows from standard arguments in kinetic theory of gases: one evaluates 

the frequency of collisions of gas particles of a given speed and subsequently integrates over all the speeds. Note that 
we have two separate contributions from the gas right (pi and Ti) and left {p2 and T2). The result reads: 



W{V\V') 



Lpi /+^ dv^iv^ - V')H [v^ - V] ^i{v,)5 
Lp2 dv^iV' - v,)H [V - v^] Mvcc)S 



v 
v 



iiV <V' 
iiV > 



with H the Heaviside function, 5 the Dirac distribution and 
T,: 



the Maxwell-Boltzmann distribution at temperature 



■ exp 



-mvt, 



2nkBTi \2kBTi 

Performing the integrals over the speed gives the following explicit result for the transition probability: 

2 



wiv\v') 



Lpi 

Lp2 



m + M 
2m 

m + M' 



2m 



{V -V)H[V' -V] (l)i{V' ■ 



{V -v')H[v -v] My' 



M 



2m 
m + M 



2m 



iV'~V)) 
{V'-V)). 



(10) 
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From the transition probability, the rescaled jump moments A„(x) ^ can be calculated. The exact expression for 
the n-th jump moment is as follows : 



(l + e2)n off 



J, n/2p 



2 + n 



P2T^ exp 



2T2 



(-l)>iTi("+'^/'exp 



2 + n 1 Toff 2 2 
■X e 



Toff 2 2 
X e 




2Ti 





2 + n 1 
2 ' 2 ' 2ri ' 



2 ' 2 ' 2T2 



23"/22, 



ri+l 



7rm(l + £2)" 

,ri/2 



y(l-n)/2p 



"3 + n' 








2 





nn/2 



2Ti 



3 + n 3 Toff 
2 '2'27\' 



x^e^ 



- P2T2 exp 
with r the Gamma function: 



^x^e^ 
2T2 



3 + n 3 Toff 2 2 



2 ' 2 ' 2T2 



1 



r[i + fc] = ki 



2fc+ 1 



2fc+i 



1-3-5' 



(2fc + l), 



and $ the Kummer function, in its integral representation given by 

T[b] 



$[a,6,z] 



T[b - a]r[a] 



6-a-l 



(11) 



(12) 



(13) 



C. Stationary speed 

The moment equations form together with the explicit expressions for the jump moments the starting point 
for a straightforward perturbation in terms of the small parameter e. To simplify notation, we introduce: 

off 

g{n) = L\\ — ^ . (15) 

Also, the limit e — )■ entails a slowing down of the motion of the piston, which can be accounted for by introducing 
a new scaled time variable: 

T = e^i. 

The equations for the first and second moment, expanded up to order and respectively, are as follows (the 
expansion for higher moments up to the sixth moment can be found in appendix A): 

dr{x) = -2 /[1/2] {x) - 9[0] {x')s + i (6 /[1/2] {x) f[-l/2] {x')) (16) 

+ g[0] {x')e' + (^-2 /[1/2] {x) + i /[-1/2] {x') + ^ /[-3/2] {x")^ - g[0] {x')e' + 0{e') 

dAx') = 4 ( /[3/2] - /[1/2] {x^)) - 2 g[Q] {x')e + 2 (^-4 /[3/2] + 5 /[1/2] {x^) - 1 /[-1/2] {x^)^ (17) 

+4 g[0] {x')s' + (^12 /[3/2] - 16 /[1/2] {x') + ^ /[-1/2] {x') + 1 /[-3/2] (a:^)^ + 0{e'). 

Note that the condition of macroscopic equilibrium {piTi = P2T2) was used to derive these equations. In particular, 
without this constraint, an additional term corresponding to a constant, velocity-independent, force acting on the 
piston, would be present in Eq. (|16|l . 

To lowest order in e, the equation for the first moment, Eq. H16|) . is not coupled to higher order moments. It 
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displays the usual linear relaxation term of the velocity, namely, in original variables, Mdt{V) 
coefficient 7: 



7 = 4L 



2ti 



92- 



-j{V), with friction 



(18) 



For Ti = T2 and (consequently) pi — P2, this result is in agreement with [1^. We conclude that at this order of the 
perturbation, the steady state speed is zero. This is not surprising since any asymmetry is buried at the level of linear 
response theory. 

Going beyond the lowest order, one enters into the domain where fluctuations and nonlinearity are intertwined. The 
first moment is now coupled to the higher order moments. Therefore, we focus on the steady state speed reached by 
the piston in the long time limit. We will omit, for simplicity of notation, a superscript st to refer to this stationary 
regime. Recalling that we defined the effective temperature Toff by the condition {x'^) = 1 at the lowest order in 
e, we immediately find from Eq. Hl()|) that at order e the piston will indeed develop a nonzero average systematic 
speed equal to eg(0)/[2/(l/2)]. The explicit value of Tcff follows from Eq. H17|l . implying at lowest order in e that 
(x^) — /(3/2)//(l/2) — 1. In original variables, cf. Eq. lO, these results read as follows: 



TeS = VT1T2 



and 



(V) 




(19) 



(20) 



Although there is no macroscopic force present (pressures on both sides of the piston are equal), the piston attains a 
stationary state with a non-zero average velocity toward the higher temperature region. Fluctuations conspire with 
the spatial asymmetry to induce a net motion in the absence of macroscopic forces. It is also clear from Ea. (|2(JII that 
the net motion vanishes when Ti = T2 and also in the macroscopic limit M — > 00. The above result was already 
derived in 14] , but the calculation presented here is streamlined so as to allow for a swift calculation of higher order 
corrections. 

At each of the next orders, a coupling arises to a next higher order moment. We shall present here the results up to 
order e^, requiring the evaluation of the moments (x^), {x^), {x"^}, {x^} and {x^), up to orders e'^, e^, e^, and e", 
respectively (cf. appendix A for details of the calculation) . The resulting expression for the average stationary speed 
in the original variable V up to fifth order in e is: 



(V) 



m \ 1/2 
M 



nks 1 
2M 2 



m) ' 



3 piVTi + P2VT2 



/T2 



(-) 



5/2 



1 



2M 

85 pi^/% + P2VT1 



V 



-1/2 



P2T.^ 



-1/2 



16 ^/T{^2 

' T1T2 {VTi ' Vn) 



+ 



18 Pl^l+P2^/% 

47. {p.T;'" + P2T2"') iVT, 

12 (pi\/rr + P2\/T2) 



Ipl^l 

3 



-3/2 



(Pl\/7Y + P2\/T2)' 

-3/2 



P2T2 



3 P1VT1+P2VT2 
3 



-T1T2 



297r (VT\ 



12 



T1T2 



(21) 



As required, the average speed is zero at equilibrium, when Ti ~ T2. Note also that the average speed depends on the 
densities of the gases solely through their ratio pxj P2- This implies that, for T\ and T2 fixed, varying the densities 
will not modify the steady-state velocity when operating at mechanical equilibrium. In Fig. 12^. and Fig. I^l^., we 
illustrate the dependence of (V) on the temperatures: the piston always moves towards the high temperature region 
and its speed increases with the temperature difference. 
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Figure 2: Stationary average speed of the adiabatic piston according to Eq. 1)21^ . Left: For pi = 0.25 and Ti — 2.0 fixed, the 
speed is shown as a function of Tb- Note that p2 is determined by the condition of mechanical equilibrium (piTi = P2T2) and 
that the piston always moves to the higher temperature region. Right: the speed of the piston as a function of Ti and T2, for 
pi — 0.25 and p2 = piTi/72- The velocity vanishes when Ti — T2 and is maximal for a large temperature difference. The 
following parameter values were used: mass of the gas particles m = 1, mass of the piston M = 100 and fc_g = 1 by choice of 
units. 

IV. THERMAL BROWNIAN MOTOR 
A. Motivation 

The systematic motion observed in the adiabatic piston is not entirely surprising since the piston is embedded in 
a nonequilibrium state with an expHcit spatial asymmetry of its surroundings. More interesting is the case of the 
thermal Brownian motor, which was introduced and studied by molecular dynamics in a recent paper [l^ . While 
the spatial environment is perfectly symmetric here, the object itself has a spatial asymmetry. The nonequilibrium 
conditions are generated by its interaction with two (or more) ideal gases that are not at the same temperature. 
The perturbative analysis, presented for the adiabatic piston, can be repeated here but is more involved because the 
problem is now genuinely two-dimensional. 



B. Presentation of the model 



Consider a 2-dimensional convex and closed object with total circumference S. Suppose that dS is a small part 
of the surface, inclined at an angle 9, measured counterclockwise from the x-axis (see Fig. O. We define the form 
factor F{0) as the fraction of the surface with orientation 0. This means that SF{9)d9 is the length of the surface 
with orientation between 9 and 9 + d9. One can immediately verify that F satisfies 



F{9) ^ 

S'^^d9F{e) = i 

J^''d9Fi9)sm9 = J^ 



positivity (a) 
normalization (b) 
d9F{9) cos 9 = object is closed. (c) 



(22) 



To simplify notation we will write (sin 9) instead of d9F{9) 

We suppose that the object, with total mass M and velocity V, has no rotational degree of freedom and a single 
translational degree of freedom. Choosing the latter oriented following the x-axis, we can write V — {V, 0). Collisions 
of the gas particles, of mass m and velocity v, with the object are supposed to be instantaneous and perfectly elastical. 
Hence, pre-collisional and post-coUisional velocities of the object, V and V, and of a gas particle, v' — {v'^,v'y) and 



(Wx 



are linked by conservation of the total energy and the momentum in the x-direction. 



1 ,2 
-mvy 



2 



1 2 1 2 



mv' + MV' = mv^ + MV. 



(23) 
(24) 



Furthermore, we assume a (short-range) central force, implying that the component of the momentum of the gas 
particle along the contact surface of the object is conserved: 



V ■ eii 



(25) 
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e + de 




Figure 3: A closed and convex object with total circumference S. The length of the surface with an orientation between 9 and 
e + de is SF{e)de, defining the form factor F{e). 



with eii = (cos 6*, sin 6'), see Fig. |21 This yields for the post-coUisional speed V: 



2^sin^ 

V = V + ^ 



1 + ^sin^ 



(26) 



Similar as in the adiabatic piston problem, we start from the linear Boltzmann equation (Q, which is exact in the 
ideal gas limit, to describe the motion of the object. The object consists out of rigidly linked (closed and convex) 
parts, each sitting in a reservoir i containing an ideal gas with uniform particle density pi and Maxwellian velocity 
distribution (pi at temperature T,;: 



2nkBT, 



exp 



(27) 



Examples of the construction with two reservoirs are schematically represented in Fig. 0] The transition probability 
M^(F|F') is then the sum of the contributions of the different units of the object and can be calculated, starting from 
the basic arguments of the kinetic theory. The contribution dWi to W{V\V') of the surface section of size dSi, with 
orientation in [9, 9 + d9], exposed to the gas mixture i is 



+ OC 



dW^{V\V') = S,F,{9)d9 



dv' 



dv'yH 



{V' ~ v') ■ e_i {V' - z/) • ej_ 



2^sin^ 
V -V — 



M 



(28) 



with H the Heaviside function, 5 the Dirac distribution and e±_ = {sin9, — cos0) a unit vector normal to the surface, 
see Fig. The total transition probability is then given by 



W{V\V') = y2 rdW,iV\V'). 
, Jo 

The integrals over the speed of the colliding gas particles can be performed explicitly, resulting in: 

( {V' - V)H[V' - V] I +{V - V')H[V - V'] 

\ Jsine>0 



(29) 



WiV\V') = 1^5, 



d9F,{9) 



2TikBTi 
M 



msvaO 



exp 



n{V' + \[{V-V') {1 + ^)]) sm' 



sine<0. 
2 



(30) 



The remaining integral depends on the shape of the object. 

Following the general set-up, we switch to the dimensionless variable x = 
temperature T^s to be determined from the condition (x^) = 1. The exact expression for the n-th rescaled jump 



M 



■V , cf. with the effective 
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(c) (d) 

Figure 4: Four different realizations of the thermal Brownian motor, each consisting of two rigidly linked units. The units are 
simple convex obj ects: a bar of length L, a disk of radius R or an isosceles triangle with apex angle 29o and base L. (a) was 
introduced in and will be referred to here as Triangula. 



moment j4„(x), cf. ITJ, is then: 



knTi ( T, 



■Km vTcff 



^/2 



-n/2 



+ (-1)"/ )d0F,(0)exp 

sin0<O Jsin0>O, 



Tcs sin 



sin ( 



r 



r n' 













n 1 TcffSin Q 2 2 



1 + 77,:^ 



2 ' 2 ' 2T, 



-X e 



ex smfi/ 1 



1 + sin^ 

+ 



3 + 71 3 Tcff sin^ 9 



2 ' 2 ' 2T, 



x^e^ 



(31) 



with r the Gamma function H12fl and $ the Kummer function (|13|l . 

C. Stationary speed 



The equations of moments together with the explicit expressions for the jump moments (|31|l . provide the 
starting point for a straightforward series expansion in terms of the small parameter e. The equations for the first 
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and second moment, expanded up to order and respectively, are given by (r = e^t): 

dr(x) ^ ' ^ . - ' f^«' 




(32) 




IksTi 




'_T}_ 


1 m 




^ TcS 


-^^ 




I TcS 
' 3 T, 


-16(x2 


7 TcS 
6 T, 


{x') + 




(33) 



Note that a term of order in the series for dr{x) is zero because of the property (1221 0). Such a term would 
correspond to a constant, velocity-independent, force acting on the object. It should indeed be zero, since each gas 
mixture separately is in equilibrium. From Eq. H32|) we also immediately recognize to lowest order in e the linear 
relaxation law, written in original variables as Mdt{V) — —^V, with 7 = 7* ^^'^ 7* linear friction coefficient, 
due to the section of the motor sitting in gas mixture i: 



7i = ASiPi 



2tt 



(34) 



At this level of the perturbation, the speed of the object is zero: no rectification takes place at the level of linear 
response. 

In order to find the first non-zero contribution to the velocity (x), we need the terms up to order s. From the 
definition of Toff, by the condition (x^) = 1 up to lowest order in e, we find from Eq. H33|l : 



TcS — 



E»7» 



(35) 



The lowest non-zero term for the average velocity then follows from Eq. (|32() and reads in the original variable V 

f^^i^EiStP, (4 - 1) !;-deF,(e)mn'e 



8M 



J:^ S^P^ J ^jTdOF.m sin' (6) 



(36) 



This speed is equal to the expansion parameter times the thermal speed of the motor and further multiplied by a 
factor that depends on the geometric properties of the object. Note that the Brownian motor ceases to function in the 
absence of a temperature difference (when Tj ~ Tcs for all i) and in the macroscopic limit M ^ 00 (since {V) ^ ^/M). 
Note also that the speed is scale-independent, i.e., independent of the actual size of the motor units: {V) is invariant 
under the rescaling Si to CSi. To isolate more clearly the effect of the asymmetry of the motor on its speed, we focus 
on the case where the units have the same shape in all compartments, i.e. Fi{6) = F{9) and Si = S. One finds: 



TcS — 



.3/2 



and 



nksTcs P^inTf- 



8M 



Ei Pi 



(sm 



(37) 



(38) 



11 



Shape Circumference S Form factor F(^) Friction coefficient 7 



Bar 2L 1 (5 _ |] + ^ _ ) 8Lp./^ 

Disk 2tvR 1/2tv AivRp^^^^ 

Triangle L1±^ 24.-^] .in.o+.i.-.oi+.i.-(.-.o)i ^rfl^^ .^^ ^ 



Table 1: The circumference S, the form factor F{0) and the friction coefficient 7 in a gas with density p and temperature T, 
for a vertical bar of length L, a disk with radius R and an isosceles triangle with base L and apex angle 29o, pointed in the 
a;-direction. 



In this case T^g is independent of F{9) and the drift velocity is proportional to (sin'^ 6')/(sin^ 6*), with the average 
defined with respect to F{d). The latter ratio is in absolute value always smaller than 1, a value that can be reached 
for "strongly" asymmetric objects as will be shown below on specific examples. The resulting speed is then very large, 
i.e., comparable to the thermal speed. 

Calculation of the average speed up to order requires the evaluation of {x'^), {x^), (2;*), (x^) and (x^), up to orders 
e^, , e^, e and e° respectively, cf. appendix B. This calculation has been carried out using symbolic manipulations, 
and the resulting expressions have been used in Table HVI and Fig. |31 However, since the analytic expressions are very 
involved, we only reproduce the result here up to order s^: 




(39) 



D. Special cases 

The above analytic result, Eq. (|39|) . is valid for any convex shape of the constituting pieces. To illustrate the 
type of explicit results that are obtained, we focus on simple shapes like a disk, a bar and a triangle. In Tabled the 
circumference S, the form factor F{6) and the friction coefficient 7 are calculated for these objects. Note that the 
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C it Lli t; 


c»^~i /~i"n c* /\/\ ( /^yri ov c i 
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Triangula 


Fig. BK. 


^27rfcBm^ ^pjp2(Ti-T2)(x/7T-v^) 
[PlV J-l+P2V-'2\ 


Triangle - bar 


r Ig. I4D. 


^2irfcBm^^ „;„2 \ 2pi pa VTr(Ti -T2 ) 

^ ''"^[2pi^+P2^(l + Bineo)]^ 


Triangle - triangle 


Fig. BJ:. 


[Pl V-'l+P2V-'2j 


Triangle - disk 


Fig. Hi. 


^27vkBrn, „;„2flX7rR PI P2 \/^(ri -T2 ) 



Table II: Analytic result for the lowest order contribution to (V^) for the different constructions depicted in Fig. 0] 



friction coefRcient of the bar is in agreement with the result of the adiabatic piston problem, cf. Eq. (|18|l . 

A thermal Brownian motor can only operate under nonequilibrium conditions, which can be achieved if at least 
two of such units are each located in a reservoir containing an ideal gas at a different temperature. The two units are 
rigidly linked and can move as a single degree of freedom along the x-direction. Besides the nonequilibrium constraint, 
a spatial asymmetry is also required to yield a net motion. In particular, a construction with only bars and/or disks 
will not generate any net motion. One of the simplest motors one can imagine was introduced in and will be 
referred to as Triangula, see figureBJi. Two identical rigidly linked isosceles triangles (with apex angle 20o, pointing in 
the x-direction) are each located in a reservoir containing a gas separately in equilibrium at a different temperature. 
The lowest order contribution to the average velocity of Triangula follows from Eq. H3t)|) , cf. Table |3 



(V^)Tda„gula = (1-Smgo) ^2 ' (40) 

The fact that the combination of the asymmetry and the temperature gradient is necessary to generate the systematic 
motion, is contained in Eq. H4()|l . If either Ti = T2 or = 7r/2 (the triangle becomes a bar and thus the asymmetry 
disappears), the average velocity vanishes. The corrections of order to this formula can be found in appendix 
B. The dependence of the speed on the temperatures and densities are reproduced in Fig. |5| Note that from Eq. 
(|40|l may be concluded that the speed of Triangula is maximal for pi^/Ti — P2^/T2- Fig. 01 also shows some other 
variations of the thermal Brownian motor. Their drift speed up to lowest order in e can be found in Table Hll One can 
easily verify that each thermal Brownian motor ceases to function when the spatial asymmetry or the temperature 
difference vanish. In this context it is sometimes stated that equilibrium is a point of flux reversal. This is indeed the 
case for our microscopic model except when there are special symmetries in the system. When the units in the two 
reservoirs are not the same, the direction of the net motion changes when the temperature difference changes sign. 
From the models depicted in Fig. 0] only Triangula keeps its original direction of motion. In this latter case, the 
speed exhibits a parabolic minimum as a function of the temperature, with zero speed at equilibrium. The reason 
for this peculiar behavior derives from the permutational symmetry of identical units, implying that the speed must 
be invariant under the interchange of Ti,pi with T2,P2- lu particular it must be an even function of Ti — T2 when 
Pi = P2- 



V. COMPARISON WITH SIMULATIONS 



The analytic results for the adiabatic piston and the thermal Brownian motor are compared with the results of the 
numerical solution of the Boltzmann equation in Table IIIII and Table IIVI respectively. To improve the precision of 
these simulations, we used a special technique for solving a Master equation, based on the introduction of a simple 
envelope process, see appendix C for details. The agreement between theory and simulations is very satisfactory, and 
only breaks down, as expected, for small ratios m/M where the perturbative result becomes inaccurate. We have 
also included for completeness the results obtained by molecular dynamics simulations, cf. T^l for more details. We 
used low densities for the hard disks, p = .0022, a regime in which one comes close to the properties of an ideal gas. 
Nevertheless one expects strong finite size effects, due to the fact that the reservoirs containing the gases are not 
very large. To cite just one important phenomenon, we note that the motion of the motor will generate sound waves 
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<v> 




Figure 5: Average speed of Triangula according to Eq. IIB4t . Left: the dependence of the velocity on the densities (Ti = 1.0 
and T2 = 5.0). Right: The stationary speed increases with the temperature difference (pi = p2 = 0.0022). The following 
parameter values were used: mass of the gas particles m — 1, total mass of the motor A4 = 100, apex angle of the triangles 
26^0 ~ n/18 and fcs = 1 by choice of units. 



Mass 


Theory 


Theory 


Theory 


Boltzmann 


M 


(order e) 


(order e + e^) 


(order e + + e^) 


equation 


1 


5.64 


-7.82 


20.44 


1.411 


5 


1.13 


0.590 


0.8158 


0.7289 


20 


0.282 


0.2484 


0.2519 


0.2511 


50 


0.113 


0.1074 


0.1076 


0.1076 


100 


0.0564 


0.05505 


0.05508 


0.0554 


200 


0.0282 


0.02786 


0.02787 


0.0280 



Table III: Stationary average speed of the adiabatic piston from the perturbative solution method up to order e, and e^, 
compared with the result from a numerical solution of the Boltzmann equation. The following parameter values were used: 
particle densities pi — 10~^ and p2 = 1, temperatures Ti = 100 and T2 = 1 and mass of the gas particles m = 1. fcs = 1 by 
choice of units. 

that will reimpact on it. Taking this into account, the speeds observed in the molecular dynamics are in reasonable 
agreement with the theoretical and numerical results obtained from the Boltzmann equation. Notably the speed of 
the thermal motor in the hard disk gases is systematically larger by roughly 20 to 40%, for reasons that are unclear 
to us. 



Mass 


Theory 


Theory 


Theory 


Boltzmann 


Molecular 


M 


(order e) 


(order e + e^) 


(order e + + e^) 


equation 


Dynamics 


1 


0.38 


-1.23 


11.66 


0.057 


0.12 


5 


0.076 


0.01188 


0.1150 


0.0470 


0.064 


20 


0.0190 


0.01502 


0.01663 


0.0157 


0.024 


50 


0.00762 


0.006974 


0.007077 


0.0071 


0.0093 


100 


0.00381 


0.003648 


0.003661 


0.0035 


0.0043 


200 


0.00190 


0.001864 


0.001866 


0.0017 


0.0021 



Table IV: Stationary average speed of the thermal Brownian motor Triangula from the perturbative solution method up to 
order e, and e^, compared with the result from a numerical solution of the Boltzmann equation. The following parameter 
values were used: particle densities pi — P2 ~ 0.00222, temperatures Ti — 1.9 and T2 = 0.1, mass of the particles m — 1 
and apex angle of the tri ang le 26^0 = 7r/18. fcs = 1 by choice of units. For comparison, we also include the speed observed in 
molecular dynamics, see jlal for more details. 
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VI. DISCUSSION 



The problem of the Maxwell demon has been haunting the imagination and theoretical efforts of physicists for more 
than hundred years. While there is a consensus that one cannot rectify fluctuations at equilibrium, it is recomforting 
that one can construct microscopic models, involving interactions with ideal gases only, for which this thesis can be 
verified explicitly. The same models can be used as test cases for another important field of interest, namely the 
rectification of thermal fluctuations in nonequilibrium, also referred to as Brownian motors. In this respect we claim 
that our model is a genuine Brownian motor: the rectification appears at the level of nonlinear response, where the 
usual separation between systematic and noise terms, as made explicit in a linear Langevin equation, is no longer 
possible. Hence the operation of our Brownian motor falls outside the scope of linear irreversible thermodynamics. It 
belongs to the realm of microscopic theory in which nonlinearity and noise form an intertwined part of the microscopic 
dynamics. This is in contrast to most of the Brownian motors discussed with mesoscopic theory. For example, the 
prototype model of one class of Brownian motors referred to as flashing ratchet can be described by diffusion in an 
external potential, a standard problem in linear irreversible thermodynamics. 



Appendix A: DETAILS ON THE CALCULATION FOR THE ADIABATIC PISTON 

Expansion of the rescaled jump moments An{x) (|11() in the moment equations ^ yields the following power series 
expansion in e = m/M, with the first and second moment expanded up to fifth and fourth order in e (cf. Eqs. 
(IIGII?!) ) respectively, and the higher moments (x^), {x"^), {x^) and (x^) up to e^,£^,e and e° respectively: 

dr{x') = 6(2/[3/2] (x)- /[1/2] (x^)) - 3 (4 g[2] + ^[O] (x^)) e 

+ (-56 /[3/2] (x) + 24 /[1/2] (x^) - /[-1/2] {x')) + 9 (4 g[2] + g[0] (x^)) + 0{e^) 
dr{x^) = 8(3/[3/2] {x?)- /[1/2] (x^)) - 4 (12 .9[2] (a;) + .9[0] {x^))e 

+4 (^16 /[5/2] - 44 /[3/2] (x?) + 11 /[1/2] (x^) - \ /[-1/2] {x^)\^ + 0{^) 

dM) = 10 (4 /[3/2] - /[1/2] (x^)) - 5 (24 g[2] {x?) + g[0] (x^)) e + ©(e^) 

dr{x^) = 12 (5 /[3/2] (x^) - /[1/2] (x6)) + 0(e). (Al) 

In the stationary regime these equations form, together with Eqs. H16I17|1 . an algebraic set of equations, which can 
be solved to find the stationary average velocity (x) up to order e^. The result in terms of the original variable 
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V — \l ^^^^x is given in Eq. H21() . The corresponding results for the higher order moments read: 
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(A2) 



Note that the lowest order terms are consistent with the observation that the velocity distribution itself is Gaussian 
at the lowest order, cf. 1141 . 



Appendix B: DETAILS ON THE CALCULATION FOR A THERMAL BROWNIAN MOTOR OF A 

GENERAL SHAPE 

A perturbative series in e for the moments, obtained by the expansion of Eqs. 0, together with the jump moments 
defined in (|11|) . The resulting expressions for the first and second moment up to order and e"* respectively are given 
in Eqs. (|32I33(I . Calculation of {V) up to requires furthermore {x^), {x'^), {x^) and (x^) up to order e'^, e^, e and 
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respectively: 
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(sm 
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(sm 
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^U^{x') + U{x')-\^{x^) 

Jeff O Ji 



(sin* 61). 



TT V Toff 



X sm 



(Bl) 



In the stationary regime, these equations form together with Eqs. I32I33|I an algebraic set of equations from which 
the average velocity can be obtained up to order The analytical expression in the original variable V up to is 
reproduced in Eq. H39|l . The corresponding power series for the higher moments of the stationary velocity distribution 
function are: 
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For the particular case of Triangula, see Fig. one can verify that 



PivTi + P2VT2 
(-1)" sin 6*0 + sin" 0o 
1 + sin do 



The average speed of Triangula up to reads: 
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m\3/2 /TrfcsTeff f 1 



2M 



sin^ 00 - 1 7r/i[0] 
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Appendix C: SIMULATION OF THE MASTER EQUATION 

Consider a Markov process, defined by the transition rate VF(T^|V^'), for transitions from state V to V. We 
decompose the transition rate as follows 

W{V\V') ^ R{V')P{V\V') 

where R{V') = J dV W{V\V') is the total rate and P{V\V') is the conditional probability. Stochastic trajectories 
of V may be easily produced if the probability P can be easily realized, that is, if it is simple to generate random 
values with that distribution. However, for complicated W and P this direct approach is often not possible and an 
alternative construction is necessary. 

We introduce an "envelope" process, W*{V\V'), for which W* > W for all Vy', cf. The difference 

Wo{V') = W*{V\V') - W{V\V') 

is called the "null" process. The envelope process is chosen such that R* is a constant and that P* has a simple form 
(e.g., a Gaussian distribution). 

The simulation of trajectories of W then proceeds as follows: 

1) Given the current state V , find the time to the next transition for the envelope process as 

T = -\n(n)/R* 

where 3? is a uniformly distributed random value in (0, 1). The random variable r is exponentially distributed with 
(r) = 1/R*. 

2) Choose the new state, V, from the distribution P*{V\V'). 

3) With probability VF(F|F')/W^*(F|y ) the transition V' ^ V occurs, otherwise, the event is a null event and the 
state V is unchanged. In either case, the time is advanced by r. 

4) Return to step 1) until the required number of iterations are performed. 

Note that the probability of a transition V' ^ V being selected in step 2) is P*{V\V') and the probability of that 
transition being accepted in step 3) is W/W* so the net probability of the accepted transition is 

'w*iV\V') 'r*P*{V\V') i?*'^*^!"^^ 
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Since the total rate (accepted plus null transitions) is R* the algorithm produces the stochastic process with the 

desired transition rate 

Clearly, the method will be inefficient if the ratio W/W* is small since most transitions would be rejected. On 
the other hand, the method is only correct if W/W* < 1 for all V, V since the probability of acceptance cannot be 
greater than one (i.e., the null process cannot have negative probability for any V). As such, the form of the envelope 
process (e.g., the mean and variance of the Gaussian) must be chosen with care. 
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